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Abstract. We analyze the diffusive motion of kink solitons governed by the thermal sine-Gordon equation. 
We analytically calculate the correlation function of the position of the kink center as well as the diffusion 
coefficient, both up to second-order in temperature. We find that the kink behavior is very similar to 
that obtained in the overdamped limit: There is a quadratic dependence on temperature in the diffusion 
coefficient that comes from the interaction among the kink and phonons, and the average value of the wave 
function increases with yjt due to the variance of the centers of individual realizations and not due to kink 
distortions. These analytical results are fully confirmed by numerical simulations. 

PACS. 05.40.-a Fluctuation phenomena, random processes, noise and Brownian motion - 05.45.-a Non- 
linear dynamics and nonlinear dynamical systems - 74.50.-fr Proximity effects, weak links, tunneling 
phenomena, and Josephson effects - 85.25.Cp Josephson devices 

rate from the zero- voltage state induced by thermal fluctu- 
ations, obtaining very satisfactory results compared with 
the experimental ones. 

Specifically, this work is devoted to the study of the dif- 
fusive dynamics of sG kink solitons subjected to a thermal 
bath, as given by the stochastically perturbed, damped sG 
equation 

(/>tt - 4>xx + sin(0) = -a^t + /(a;, t, (/),...), (1) 

where —a(f)t is a damping term with a dissipation coefR- 
cient a, and f(x, t, (j)^...) is a thermal (gaussian) noise term 
fulfilling 

f{x,t,(j)^...) = VD r]{x,t), {r]{x,t))=0, 

(2) 

{rtix, t)ri{x',t')) = S{x - x')Sit - t'), 

where ^/D is related to temperature through the fluctu- 
ation-dissipation theorem D = 2akbT, kf, being the Boltz- 
mann constant and T the temperature. 

To our knowledge, the flrst results on problems di- 
rectly related to the one we deal with here were obtained 
by Joergensen et al. |26|| , who performed experiments on 
Josephson junctions and presented a derivation of the dif- 
fusion constant for kinks. Subsequently, Kaup and Os- 
man | |2^ ] studied, in a more rigorous way, the motion of 
damped sG kinks, driven by a constant force, in the pres- 
ence of thermal fluctuations by using a singular pertur- 
bation expansion. They analyzed the temperature effect 



1 Introduction 

As a key subject within nonlinear science, the dynam- 

■ ics of emergent, coherent structures (solitons, vortices, 
\ etc) has been a research topic that has attracted very 

■ much attention in the past quarter century [Q . One ques- 

■ tion, extensively investigated in the literature |^,|^,^,||,||, 
1^ is the following: Is, and if so, how is the motion and 

' the shape of those excitations modified by the presence 
' of small perturbations? Indeed, when applied to physical 
situations of interest, nonlinear models must incorporate 
' additional terms, such as damping, constant or periodic 
, external forces, or noise, to name a few. Among those, 

■ stochastic perturbations are very much of interest in view 
' of their highly non trivial effects on nonlinear system s H , 
. and a great deal of work has been devoted to them 

In particular, of the very many nonlinear models ap- 
, plied to physical problems, the sine-Gordon (sG) equation 
has been considered in much detail in this context, as it 
applies to, e.g., propagation of ultra-short optical pulses 
in resonant laser media 0], a unitary theory of elemen- 
tary particles ]lO|, pd] , [l2i|l^ , propagation of magnetic fiux 
in Jose phs on junctions ||14||, transmission of ferro mag netic 
waves epitaxial growth of thin films ||l^,|l^,|8||, mo- 
tion of dislocations in crysta ls [p^[20[| , flux-line unlocking 
in type II superconductors |21|, or DNA dynamics p2| , 
situations in which noise (of different origins) can 
play, and often does, a crucial role. As an example, let 
us mention the recent work on long Josephson junctions 
reported in |23], where the authors calculated the escape 
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on the mean velocity of the kink and also the changes 
in the shape of the kink. In addition, they calculated the 
diffusion coefhcient of the kink up to first-order in temper- 
ature and the energy values corresponding to the transla- 
tional {Et = ki,T/2) and radiational {Eji = kbT) modes. 
These values of the energy have been also obtained by 
Marchesonip^ , who applied the McLaughlin and Scott 
approach [|29|| to investigate kink motion under thermal 
fluctuations (see |||,||,|| for reviews). 

For the sake of completeness, let us mention work done 
along a different line, namely that devoted to the diffusive 
motion of the kink in equilibrium with phonons in isolated 
sG systems (possibly perturbed) |||,^,^,32 . In this case, 
the kink diffusive motion is characterized by two diffusion 
coefficients. The first one of them is proportional to 
and is related to the anomalous diffusion, that arises from 
the phase shifts of kinks colliding with phonons and takes 
place on a short time scale in which the collision among 
kink and wave packet is elastic; the kink retains the same 
velocity after the collision (non-dissipative diffusion) and 
suffers only a spatial shift. However, for large times and in 
slightly perturbed sG systems, this interaction is nonlinear 
and becomes inelastic and the velocity of the kink changes 
after the collisions |Q . This diffusive regime is called vis- 
cous and has a diffusion coefficient proportional to T^^. 
The diffusion of the kink when the low energy excitations 
are rep resented by breathers has also been studied, and 
in |34| it has been demonstrated that both descriptions 
(breathers or phonons) are equivalent and give rise to the 
same diffusion coefficient in the anomalous regime. 

In any event, we want to stress that, although in this 
type of diffusion problem there are many open questions 
p3{ , we will concern ourselves with the other kind of dif- 
fusion problem, in which the phonons appear as a con- 
sequence of an external heat bath, represented by white 
noise correlated in space and in time and the damping is 
included explicitly d la Langevin. The main aim of this 
work is to extend a previous study of ours about the over- 
damped limit of sine-Gordon kink diffusion |^ to the 
more physical and general case of the underdamped dy- 
namics (i.e., with finite dissipation coefficient). As we will 
see below, the general perturbative approach p3[ we re- 
sorted to in the overdamped case can also be applied, al- 
beit with more difficulties, to the underdamped problem. 
The corresponding theoretical analysis is presented in Sec. 
2, where we obtain explicit expressions for the long-time 
diffusive dynamics of kinks up to second-order in temper- 
ature, thus going beyond the currently available knowl- 
edge. The accuracy and importance of the new terms is as- 
sessed by numerical simulations in Sec. 3: we will see there 
that the quadratic corrections are in good agreement with 
the simulations and, most importantly, that they must be 
taken into account even for not so large temperatures. Fi- 
nally, in the conclusions we summarize the main results of 
this work, comparing the underdamped and overdamped 
dynamics of the sG equation and discussing other related 
questions. 



2 Analytical results 

We begin by briefly reviewing the basic results we need for 
our analytical approach. We will concern ourselves with 
the perturbation effect on the kink solutions of the unper- 
turbed sG equation, whose static form is 



ho{x,t) — 4 arctan[exp(j:)]. 



(3) 



Small perturbations over this equation can be treated by 
calculating the spectrum of linear excitations around the 
kink solution |3q| : To this end, we write 

(t)ix,t) ^ (j)oix) +'tP{x,t), 'ilj{x,t) «Mx), (4) 

substitute in (|l]) (with a = D = 0) and linearize around 
4>o{x), arriving at the following equation for il){x,t): 



iptt = ^xx - 



1 - 



cosh (x) 



(5) 



Then, assuming that the solution of (||) has the form 

iplxjt) — f{x) exp (^iut^ (6) 

we find the eigenvalue problem for f{x), 

2 



dx^ 



1 - 



cosh (x) 



(7) 



This equation admits the following cigcnfunctions with 
their respective eigenvalues 



cosh(a;) 



LU^ = 0, 



(8) 



fkix) = 



exp^ikx) [k + itanh(a;)] 



2tt iOk 



, ujl = l + k\ (9) 



which represent, respectively, the translation (Goldstone) 
mode and the radiation modes. Importantly, the functions 
frix) and fk{x) form a complete set with the orthogonal- 
ity relations 



+ OC 



fx{x) dx 



+ 00 



fT{x)fk{x) dx = 0, 
fk{x)f;,{x) dx = S{k-k'). 



(10) 

(11) 
(12) 



We can now proceed with our problem: In order to 
tackle Eq. (|l|) , with noise as given in (||) , we use the same 
Ansatz proposed for the overdamped case in [^5| (or for 
the general Klein- Gordon system in jl^): We assume that 
the solution of Eq. (|l|) is 



dkAk{t)fk[x-X{t)], (13) 



(x,t) = </)o[x-X(t)] + 



where X{t) is the kink position. We now insert ( |13[ ) in 
(P and use the orthogonality of fk and fr [0 , obtaining 
the following system of differential equation for X{t) and 
Ak{t): 
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dk / dk'Ak{t)Ak'{t)R^{k,k') 



/-\-oo -1 
1 
dkAk{t)h{k) - — 
-OO ^ ^ —OQ ^ — OO 

+ hv-x(t)\^{x,t) dx- 

x: /'+C)C /'+00 

dk / dfci / dfc2Afe(i)Afe,(t)Afe,(i)i?6(fc,fcl,fc2) 

'— OO J — OO J —OO 



1 

48 



dkAk{t)h{k) 



(14) 



where 



+ 00 

+ c^2^fe(t) = / dk'Ak'{t)h(k',k) + 



dt^ dt 

1 

+ 2 



r'+oc 

dfc / dk'Ak{t)Ak'{t)R4{k,k') 



D mx-X{t)] r^ix,t) dx + 



— OO 



dfc' / dfci / dk2Ak'{t)AkAt)Ak2{t)R7{k',k,ki,k2) + 



1 

/ + CXO Q A 

dk'^h{k',k) + 

r + co 

Xit) dk'Ak'{t)h{k\k) + X\t)h{k), (15) 



" A/27r cosh 

4ci;fcWfc' sinh f 



Ri{k,k')^ / [/*(0)]2_^/,(0)d0, i?4(fc,A:) 



' 8V2^ cosh(!^ 

/ + 00 
cos(0o)/fe' {0)fk{e)fl {e)h, {9)d9. (16) 
-C30 

It goes without saying that these equations can not be solved. Therefore, in order to extract information from 
them, we resort to a perturbative approach assuming the noise term is small, or equivalently, that the temperature 
and the dissipation arc not too large (this is not a serious restriction since our single-kink approach does not apply to 
high temperatures, when kink-antikink pairs are thermally generated |37| ). We then expand X{t) and Ak{t) in powers 
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of /D, i.e., X{t) = J2'^^i{VDTXnit) and Ak{t) = J2^^j^iVDTA'^{t), since when /D = and a = we recover the 
static kink sohition (in this case initially centered at the origin) of the sG equation. By substituting these expansions 
in (|lj) and (|l5|) we find a set of linear equations for the coefficients of these series. The first members of this hierarchy 
correspond to order 0{^/D): 



Xi{t) + aXi{t) = 



+ OC 



fT[x~X{t)] f]{x,t) dx = ei{t), 



(17) 



from where we obtain the statistical properties of ei{t), 



(18) 



and 



+ 00 



nix-Xit)] 7?(x,t) dx = Ut), 



(19) 



which in turn leads to 



(20) 



Equations (|l7|)-(p0|) have been obtained in |Q by using a similar, but more restrictive perturbative approach. By 
integrating these two equations we obtain the first-order terms, Xi{t) and A^: 



Xi{t) = I e 
'o 



C\it) 



-at' 



e"^ ei{T)dT dt' , Al{t)^e'2 {Ci{t)iiinujt + C2{t) cosut} , 



(21) 



S,k{T)e 2 cos o^rdr, C2{t) = / $A:(T)e 2 sincjrdr. 



where = — (a^/4). From these relations we can calculate the mean values and correlation functions up to first 
order in ^/D: 



{X,{t)) = 0, {X{t)X(t')) = D{X,{t)X,{t')) = 



D 



16a3 



-a\At\ 



-"(«+«') + e-"t + e-"t' + 2{aM - 1) 



(22) 



{Mt)) = 0, {mX{t')) = D{X,{t)X,{t')) = 

16a 



f,-a\At\ _ g-Q(t-Ht') 



(23) 



{Alit)) = 0, {A,(t)A,{t')) = D{AUt)AUt')) 



D 



'a{t+t')/2 



ae 



aM 



8-^ 



aM 



■ sa\Lij\At\ 



2a 
a 
8w 



■ coBLoAt 



2 coscj(i + 1') 



- — 7j smuit + t ) 
4wf 



(24) 



where At — t~ t', and M = min(t, t'). Of course, for t' = t in Eq. ( p2| ) we recover the result in [g^ for {[X{t)]'^). 

We now turn to the main point of our work, namely obtaining the next-order corrections for the position and the 
velocity of the center of the kink. This requires the calculation of the next two contributions to X{t) as well as the 
second order in the radiation terms, which are: 



0{D) 



X2(t) + aX2{t) = e2{t), 



(25) 
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dk' / dkAl{t)Al,{t)R3{k,k'), (26) 



8 .L^. «w ' 4 ./_^ dt 

1 

16 



-^{t)+a-^{t)+ulAl{t)^e,{t) / dk'Al,{t)h{k',k) 



0{WD] 



1 p + co p + oc 

+ -/ dk dk'Al{t)Al,{t)R^{k,k')- 

/ + 00 OjAl 
dk'^{t)h{k',k); (27) 



X3{t)+aX3{t) =^3{t), (28) 



£3(0 ^ r^dkAlmik) - ^ / &Ai(i)/i(fc) 



^ ^ — OO ^ 

-y dfcy dfc'A^(t)Ai,(t)i?3(fc,fc')- 

-/ dfc / dk'A\(t)Al{t)R:,{k,]^)- 

J —OO J — OO 

-/ dk dk^ dk,Al{t)Al{t)Al{t)Re{k,k,,k2) + 

J —QO J — OO J — OO 



dt ' 4 9i 

^2(0 dk^(^t)hik)- (29) 



9i 

Analogously to what we have done for the first-order corrections, from the solutions of Eqs. (25) and (Eq) we find 



that 

{X2{t))^0, {X2{t))=0, (30) 
(X3(i))=0, {X3{t))=0- (31) 

As for higher moments, taking into account that the cross-correlation function of Xi{t) and X^it') is of the same order 
as {X2{t)X2{t')), and also that {Xi{t)X2{t')) = we obtain that 

{[x{t)r) = D{[x,{t)r) + 

+ D^{{[X2{t)r)+2{X,{t)X,m + ..., (32) 

{[x{t)r)^D{[Mt)r) + 

+ D^iiiMt)?) + 2{X,{t)X,m + - (33) 

The expressions for the functions {[X2(t)]'^), {[X2{t)]'^), {Xi{t)X3{t)), and {Xi(t)X3{t)) can be obtained after a 
lengthy calculation, and are very cumbersome indeed. We therefore do not include them here. However, for large time 
{t >> 1/a) these relations can be simplified, yielding, as t — + oo 
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khT , 



9^2 
4 



[X{t)f 



with 



+ OC 



dk 



LUk cosh (-^) 



= 1.62386. 



(36) 



To complete the characterization of the kink diffusion, 
we can now compute in a straightforward way the average 
value of the wave function (f){x,t), defined as 



t)) = (M^ ~ VDXi{t)]) + 0{D) 



+ 00 



dX^p{X^)Mx - VZ?Xi(t)], (37) 



where p{Xi) is the probabihty distribution function for 
Xi. To find exphcitly this function we note that, if we 
rewrite Eq. ( |l7|) as a system of two differential equations, 



^1 
V 



-aV + ei{t), 



(38) 



the last equation represents an Ornstein-Uhlenbeck pro- 
cess for the velocity, and its distribution function is given 

by 



p{V) = 



1 



27r(F2 



■ exp 



2(1/2 



(39) 




Subsequently, by integrating the first equation 
I), we obtain that Xi = V{T)dT. Since V has a 
Gaussian distribution function, Xi has also a Gaussian 
distribution function, given by (recall that {Xi{t)) = 0) 



P(^i) = 



1 



■ exp 



(4 



X! 



(40) 



where the first and second moments of Xi were obtained 
before, see Eq. (^). With this result, the integral ( ^7|) , 
can be evaluated numerically taking into account Eqs. ( pt)|) 
and (p2[). In the next section we will compare this result 
with the mean value of the wave function as obtained from 
simulations of the full partial differential equation m). 



3 Numerical simulations 

In order to test the approximate theory developed in the 
previous section, we have simulated numerically Eq. (^ 



by using the Heun method pO| . In our simulations we be- 
gin with a kink, initially at rest, with free boundary con- 
ditions. For the damping coefficient we choose a — 0.1, 
which is not too small because from (jsj) we can see that 
([X(i)]2) is proportional to 1/a. This means that if a is 
too small the kink can move in a much larger region, forc- 
ing us to increase the length of our simulated system in the 
simulations, already quite time consuming. Furthermore, 
the relation (34) is only valid for large times (t >> I /a). 
Again, for too small a we would need to simulate our equa- 
tion for very long times and, as ([X(t)]2) increases linearly 
with time [see Eq. (|3^)], the system length would once 
more have to be large. The other parameters are Ax = 0.2, 
At = 0.001 and the length of the system L = 200. We have 
calculated all average values over 1000 realizations up to 
a final time 400. It is important to point out that, this 
system being inertial, the accuracy of the averages is con- 
siderably less than for overdamped problems, this being 
the reason why we have to use such large ensembles of 
trajectories to obtain reasonably good results. 

An important, nontrivial issue is the question as to how 
can we find the center of the kink. We solve this problem 
by finding all the discrete lattice points Xi and Xi+i such 
that 0i < TT and > t or vice versa, and then esti- 
mating the corresponding points Xi where = tt by linear 
interpolation. Afterwards, among the n such points i„, we 
choose to be the center of the kink the value x = £„, which 

minimizes '^fj,^^[4>i{t) — (j^oix — Xn, t)]'^, i.e., the discrete 
version of the integral of the square of the difference of 
(p and (j)Q. It has to be realized that this involves an as- 
sumption, namely that individual realizations of the kink 
have a shape similar to that of the unperturbed kink. As 
can be seen from Fig. ^ where the individual realizations 
are compared with the initial condition (an exact kink), 
this is indeed the case and our procedure is truly sensi- 
ble. Therefore, we are sure that this method to compute 
the kink center avoids the false centers, which can appear 
for higher temperatures due to fluctuations introducing a 
systematic difference between the numerical and the theo- 
retical results (see ^^). With the procedure we have just 
summarized, that works even for relatively large temper- 
atures, we believe we find a very accurate approximation 
to the actual center of the kink. We will come back to Fig. 
|l| below. 

As an example of the comparison of the numerical sim- 
ulations of Eq. (^ with the theoretical results obtained 
in the previous section and valid for large times. Fig. |^ 
shows the numerically computed variance of the center 
of the kink, {[X{t)f) - {X{t)f, as well as the first- and 
second-order analytical expressions. The plot clearly ev- 
idences that the numerical variance asymptotically coin- 
cides with the second-order expression: Note that to com- 
pare the different curves one has to look at the slopes at 
times t >> 1/a (in this case, t > 100, for instance, as 
a = 0.1); the theoretical result is not valid at early times 
and therefore there is a bias between analytics and numer- 
ics coming from that. The small, irregular oscillations in 
the numerical curve arise from the difficulty in accurately 
computing averages in an underdamped system like this 
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Fig. 1. Individual kink realizations compared with initial con- 
ditions and averages. Shown are three individual realizations 
a.t t — 300 with parameters ksT = 0.4 and a — 0.1 (thin solid 
lines), the initial condition given by a kink at rest (dashed 
line) and the mean value of (j) at the same time obtained from 
averaging over 1000 realizations (thick solid line). 



Fig. 2. Comparison of analytical and numerical results for the 
variance of the kink center, {[X{t)f) ~ {X{t))'^, for kt,T = 0.4. 
The dot-dashed line is the result of the numerical simulation of 
Eq. (|l]), whereas the dotted and the solid lines are the analyt- 
ical results for the first- (dotted line) and second-order (solid 
line) expressions [Eqs. ( p^ ) and (^) with t = t' respectively]. 
Only the slopes of the lines for t >> 10 have to be compared. 



mentioned above; however, we believe that the present ac- 
curacy is enough to confirm the validity of our approach. 
We have observed the same agreement for other values of 
temperatures (fcfcT = 0.2, 0.6, 0.8, not shown). In all cases, 
we have computed the diffusion coefficient for large times 
as the slope of the variance of X{t) again for t > 100, 
the regime in which we expect our analytical approxima- 
tion to be valid. Summarizing our results, these numerical 
values of the diffusion coefficient are plotted in Fig. || to- 
gether with the theoretical results. It is clear that for large 
temperatures the quadratic behavior in fcfcT of the diffu- 
sion coefficient becomes important. For higher values of 
the temperature, such as k^T = 0.8, the numerical value 
of the diffusion coefficient is not so close to the predicted 
one. This effect arises because of the large diffusivity of 
the kink in that range: Indeed, for this and higher tem- 
peratures the kink performs very long excursions away 
from the center, reaching the boundaries of the numeri- 
cal integration interval; it is clear that when this occurs, 
the diffusion of the kink is not in free space anymore and 
hence those realizations spoil the quality of the averages. 
The way to solve this problem would be to resort to much 
larger numerical systems, but within our present comput- 
ing capabilities this would necessitate a simultaneous de- 
crease in the number of realizations in the average, leading 
again to poorer results. However, it is important to real- 
ize that this boundary effect leads to an underestimation 
of the diffusion coefficient (as the boundary prevents the 
kink from travelling as far as it should) and therefore the 
point in Fig. | for fc^T = 0.8 is a lower bound for the dif- 
fusion coefhcient, with the actual one lying even closer to 
our second order prediction. 



Finally, there is one last question that deserves discus- 
sion, namely that of the physical significance of the mean 
value of the wave function 0. In Fig. |l| we can clearly see 
that, whereas individual realizations of kinks look very 
similar to the unperturbed ones, the mean value of (j) 
is a much wider excitation, not even close to the origi- 
nal kink. Figure ^ clearly shows that this does not mean 
that the width of individual kinks increases; indeed, much 
as we discussed regarding the overdamped problem [^ , 
we have verified numerically that the mean wave func- 
tion {(j){x,tfix)) increases due to the variance of the kink 
center of individual realizations, and hence it should not 
be interpreted as the typical deformation of the shape of 
kinks: Indeed, the widening of the mean value of (j) arises 
from the contributions of the stochastically moving, but 
mostly undistortcd kinks whose center positions have the 
distribution of a rigid, diffusing particle. To further check 
this interpretation, we can look at Fig. ^, where we have 
represented the mean value of the wave function for two 
fixed times tjix = 100,300, obtained from the numerical 
simulation of the full partial differential equation (|l|), for 
fcfjT — 0.6 and a — 0.1. The ovcrimposed points, com- 
puted by using the Gaussian distribution function p{Xi) 
[Eq. (|^)] of the kink center X{t) = ^/~DXi found in the 
last section, show the excellent agreement between our 
theory and the simulation. Of course, there is a small dis- 
crepancy that is likely to disappear if one would go to 
a next order calculation, but for the present purpose of 
understanding the mean wave function (j) the first order 
calculation is enough. In addition, we have plotted the 
initial kink (at rest) in order to see that the mean value 
of the wave function increases with time. 
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Fig. 3. Kink diffusion constant vs temperature. The solid lines 
represent the analytical diffusion coefficient up to first- (lower 
line) and second-order (upper line). Diamonds stand for nu- 
merical values of the kink center diffusion coefficient, obtained 
by numerical integration of Eq. (|l|. 
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Fig. 4. Mean value of the wave function (solid line) for two 
fixed times, tfi^ = 100 and tfi^ = 300, with kbT = 0.6, ob- 
tained from numerical simulations of Eq. (|l]) compared to the 
values of ((l}(x. t)) (points) obtained numerically from the in- 
tegral in Eq. (pTI). The narrowest solid line is the initial data 
(kink initially at rest). 



4 Discussion and conclusions 

To summarize, in this work we have studied the diffu- 
sive dynamics of sine-Gordon kinks subjected to thermal 
fluctuations. We have analytically calculated expressions, 
valid up to second order in temperature, for the average 
position and variance of the kink center, as well as for the 
mean shape of the kink. We have numerically checked the 
validity of these results up to temperatures of the order 



of kbT ~ 0.8 (in dimcnsionlcss units, equivalent to about 
a 10% of the kink rest mass), already close to the temper- 
ature at which kink-antikink nucleation becomes a likely 
event. Therefore, our first conclusion is that the second- 
order theory developed here is the proper one, meaning it 
is accurate and higher order terms are negligible, to de- 
scribe the thermal diffusion of sine-Gordon kinks in the 
single kink propagation regime. Interestingly, our calcu- 
lation pinpoints the fact that the second-order correc- 
tion in khT comes from the interaction between kink and 
phonons. This implies that the physics behind this contri- 
bution is basically the same as for the case of anomalous 
diffusion in an isolated chain mentioned in the introduc- 
tion |^,|o[|3l|,^. Note that we do not expect contri- 
butions in our analytical calculations, as they are carried 
out in a continuum sG equation |^ and, in any case, they 
would show up in simulations only for very low tempera- 
tures. Apart from that, it is also interesting to note that, 
according to Eq. (35), the second order term implies an in- 
crease of the energy carried by the kink beyond the k^T /2 
predicted by statistical mechanics (recall that the kink 
mass is 8 in our units). This can be interpreted in the fol- 
lowing way: The kink is dressed by phonons which increase 
its mass. Thus, the kink energy would be M{T) ■ 
with a temperature dependent mass M{T) whose expres- 
sion can be easily found from Eq. (35). In order to confirm 



this interpretation, one could compute the energy carried 
by the phonons which dress the kink, but we believe that 
it is not necessary because, on the one hand, it would be 
a rather involved calculation (far beyond the scope of this 
work) and, on the other hand, we do not think that there 
is any other possible interpretation of this result. 

A second relevant point of this study relates to the nu- 
merical simulation and center location procedures. As this 
is an underdamped (inertial) system, the thermal mobility 
of kinks is quite large, the larger the higher the tempera- 
ture. Because of this, we have not been able to obtain very 
precise numerical averages at the top of the temperature 
range studied, since the lengths of the systems and the 
number of realizations required are very large and con- 
sequently time consuming. However, we believe that the 
results presented here are enough to verify our theory. 
This is reinforced by the very good agreement between 
analytics and numerics regarding the mean shape of the 
field, even for temperatures as large as fc^T = 0.6 (see 
Fig. ^), which shows that our approach indeed captures 
the physics of the diffusion process. In addition, we want 
to emphasize that, to our knowledge, we have designed a 
new algorithm to detect the kink center which gives very 
good results even for the highest temperatures studied, 
where previous researches, such as [^^, had found prob- 
lems arising from the many false centers detected. 

Another important issue is the comparison of the pre- 
sent analysis to that in for the overdamped prob- 
lem. We have found that the diffusion coefficient given by 
( p4[ ) for the present case practically coincides with that 
obtained in |3^ for the overdamped limit of the equa- 
tion: the difference in the second order is approximately 
0.06fcf,T, i.e., very small compared to the magnitude of 
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the quadratic contribution itself. Furthermore, the width 
of the mean value of the wave function increases with 
time for the overdamped case in the same manner 
as that reported here. Therefore, we can conclude that for 
large times the dynamic of the underdamped sG kinks is 
very similar to the overdamped case. This is an important 
point, because in principle one can expect similar results 
for other kink-bearing systems such as the 0"* equation, for 
instance, whose overdamped diffusive dynamics is known 
(see for the case), thus avoiding the much more 
involved calculation of the underdamped case. 

Finally, we want to mention the relevance of this work 
to experimental systems, such as long Josephson junc- 
tions. As has been shown in the thermal sG equa- 
tion (|l|) is a good description of the physics of in-line 
Josephson junctions (although different boundary condi- 
tions are needed in that case). The work in compared 
the predictions from the sG model to experimentally mea- 
sured escape rates from the zero voltage state. There- 
fore, it should be possible to design similar experiments 
in order to test our results and, specifically, the increased 
(quadratic) diffusivity of kinks at higher temperatures vs 
the linear behavior at lower ones. We hope that our the- 
oretical work stimulates further experimental research in 
that direction. 
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